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The objective of this study was to analyze the heat transfer and fluid flow within 
a closed cavity, which usually enhance by application of Nanofluids. The move- 
ment of the Nanofluids within the enclosure were governed by Buoyancy force 
and Gravitational force. Different kinds of Nanofluids were employed for pur- 
pose of examining the motion of heat flow and the stream line formation, within 
the close cavity. The closed cavity considered for this study comprised of a flush 
mounted heater on the left wall, the right wall was maintained at a constant 
temperature throughout the investigation and the top & bottom wall were kept 
insulated. In this study, MATLAB was used for performing the simulation and 
solving the fundamentals of Computational Fluid Dynamics (CFD), namely 
the Energy equation and the Momentum equation. The analysis to estimate the 
heat transfer and fluid flow in a closed cavity by the use of different Nanofluids 
have been presented in terms of several parameters such as Rayleigh number, 


Nusselt number, Aspect ratio and Volume fraction respectively. 


1. Introduction 


Heat transfer through a natural convection mode in 
an enclosure had been in research from the past 
decades and the reason behind where the natural 
convection is the frequent process in nature itself, it 
has many industrial applications too (Azad, Groulx, 
and Donaldson Moria Kumar, Chanakya, and Bart- 
wal). During this time period, nanotechnology gain- 
ing the interest of the researcher more. Study of fluid 
flow and heat transfer by use of nanofluid helped in 
enhancing the result and value of work was compar- 
atively more efficient, presently. 


In this paper, a finite element simulation has been 
performed to inspect natural convection heat trans- 
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fer in a moderately heated rectangular cavity filled 
with nanoscale particles and base fluid as water with 
uniform thermal boundary condition (Taghilou and 
Khavasi). The enclosure was completely insulated 
at top, bottom and left wall, and the left wall was 
climb up by a flush mounted heater at the finite 
extent (Ho, Jang, and Lai). Heat transfer and fluid 
flow inside the enclosure was due to buoyancy force. 
Different volume fraction, Rayleigh number, mean 
Nusselt number and local Nusselt number was cal- 
culated for different nanofluids Cu, Ag, Pure water, 
Ti02, Al20O3. The increase in average Nusselt num- 
ber value for the change of nanofluid indicated the 
increase in heat transfer (Okabe et al. Susantez, 
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Figure 1 (a): Geometry enclosure 
mounted with the heater 
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Figure 1 (b): Velocity Boundary Condition 
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Figure 1 (c): Temperature at Boundaries 





FIGURE 1. Representing theGeometry and 
Boundary Conditions (a) Geometry enclosure 
mouthed with the heater (b) Velocity Boundary 
Condition (c) Temperature at Boundaries. 


Caldeira, and Loiola Mahdi et al.). 

Study was done by solving the non-dimensional 
Navier stokes equation with the assumption of 
Boussinesq model and finite volume method was 
used (Venkateshwar et al.). The energy equation 
was in expression of temperature, and pressure as 
the source term which was solved with different 
approach of CFD to analysis the flow phenomena 
streamlines and temperature profile by simulation 
on software, MATLAB (Ahn et al. D. Yang et al. 
Loenko, Shenoy, and Sheremet). This study talks 
about changed in heat transfer rate with respect to 
different nanofluid and with different volume frac- 
tion of nanofluid (Ma et al.). This stimulation work 
with the help of some particular solving approach of 
CFD figure out the temperature change of nanofluid 
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near the heater and as the particle near to heater get 
heated up early its molecular momentum increases 
and the particles start moving due to buoyancy 
force, the upward thrust and further new particles 
was forced to come near the left wall mounted by 
heater and heated again by same phenomena and 
so on (P. Yang et al. Khan et al. Pranowo and 
Wiyayanta) .This way the particles able to move for- 
ward toward the top wall and under the influence 
of gravitational force the particles move downward 
after reaching the end of top wall, this clockwise 
circular streamline path was seem in the particular 
enclosure (Moria). 


The primary aim of this research study was 
to solve the non-dimensional governing equation 
(Boussineq model) of CFD for natural convection 
heat transfer in a partially heated rectangular enclo- 
sure filled with nanofluids (Barman and Rao). The 
main literature review discussed further specified 
that there was no study on natural convection of 
a rectangular enclosure, for solving non dimen- 
sional governing equation of CFD by simulation 
method (by software work as a mechanical engi- 
neer) (Akinshilo Kumar, Chanakya, and Bartwal 
Alizadeh et al.). 


2. Methodology 


The Figure-l(a). shows the Geometry enclosure 
mounted with the heater. In this particular enclo- 
sure all the pressure, temperature, velocities etc was 
solved. Buoyancy force was responsible for fluid 
flow as the initial movement of particle was gen- 
erated because of the upward thrust applied to the 
heated particle near to the heater and secondly, the 
motion of fluid in the cavity is because of increase 
of temperature as temperature increases the momen- 
tum increases and fluid particles start moving and 
here the change of density of fluid was seen with the 
increase of temperature which was assumed as con- 
stant because of the approach of Boussineq model. 
We were using dimensionless study to simulate the 
work. Maximum temperature can be seen near to 
the heater so maximum velocity of fluid particle was 
near to the heater. 


The gap in research study that mainly identified 
was: 

e The numerical and experimental analysis were 
done, but the simulation work was not performed by 
any researcher. 
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e The study was not solved with the help of Com- 
putational Fluid Dynamics (CFD) approach. 

e The study was always solved under the assump- 
tion that pressure change was constant, but use of 
CFD approach to solve the pressure changes was 
never employed. 


2.1. Boundary Conditions: 


On left wall (heater): u=0, v=0, T=TyOn bottom 
wall and top wall: u=0,v=0,y=0 


2.2. Governing Equation 
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2.3. Numerical Method 
2.3.1. Finite Volume Method 


The solving approach was done by using finite vol- 
ume method of CFD to discretise the equations. We 
have three regions here: 

1- above the heater, 2- domain of heater, 3- below 
the heater 

We using Lyl, nyl for region |. Ny1 represent all 
mesh in y direction for region 1, and respectively. 








=, 1 =e Ly3 = 
= — nx—1 _ dy3 7” mya CY _ 
ei dy2 _ al 


2.3.2. Formation for P, U, V, T 


Figure-2(b). illustrates General Mesh, if three 
regions are combined then new mesh Is created and 
we reach this mesh. This mesh is very less and only 
for better understanding. We divided y direction to 
ny= 5 and x direction to nx=5. The nx and ny is 
too important because they help us for numbering of 
u velocity element and v, P, T. In the next slide we 
Show all element for u, v, P, T. 

In Figure-2(c). Mesh for U Velocity is shown 
and black vectors are boundary element and they 
are known. But blue vectors have been calculated. 
U matrix is (ny+1, nx) matrix (blue arrow position). 
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Figure 2 (a): Mesh Creation 
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Figure 2 (e): Mesh for Pressure, p 


Figure 2 (b): General Mesh 
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1.23.45 Figure 2 (f): Mesh for Temperature, T 


FIGURE 2. Mesh for Temperature (T), Velocity 
(V), Pressure (P) (a) Mesh Creation (b) General 
Mesh (c) Mesh for U Velocity (d) Mesh for Veloc- 
ity (e) Mesh for Pressure, p and (f) Mesh for Tem- 
perature, T. 





FIGURE 3. General Mesh for U, V, T, P 


Black vectors are boundary element and they are 
known. But blue vectors must be calculated. V 
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matrix is (ny, nx+1) matrix. 

Black vectors are boundary element and they are 
known. But blue vectors must be calculated. P 
matrix is (ny+1, nx+1) matrix. We add more node 
for pressure for making essay programming. All 
black node pressure is zero. Black vectors are 
boundary element and they are known. But blue 
vectors must be calculated. T matrix is (ny,nx) 
matrix. 


2.4. Boussinesq Model: 


To solve the natural-convection flows, Boussinesq 
model was used in this problem as according to this 
model the density was the function of temperature 
and the source term used is: 


(9 — po)g © — poB (T —TO0)g 


2.5. Momentum Equation: 


Pnf 


epaBs (1-9) Ps Bf )9 


dv dv _ _OP 4 Mnf ( O?u 4 Ov 
Lays + Oy Oy (e: alr ad 
+ Pnf 


2.6. TDMA Method: 


The Technique for rapidly solving tri-diagonal sys- 
tems that is now called the Thomas Algorithm or 
the tri-diagonal matrix algorithm (TDMA). TDMA 
algorithm was a two-step method that have features 
both of direct and iterative methods. iterative meth- 
ods such as Gauss-Sidle and Jacobi methods and 
direct methods such as Gaussian elimination and LU 
methods (Chen et al.). In this method, first, all vari- 
ables assumed to be unknown in one direction and 
known in the other. The second stage was the oppo- 
site of the previous stage. In this study, first the 
variables in the one direction were assumed to be 
unknown, so for the first step we have: 


1 1 1 
stepl : aw? + apU, ? = acU,,? 


= asUg + anU;, + su 


k+1 k+1 k+1 
—asUc" + apUp" — anu 


k+4 k+5 


step2 : 


2021, Vol. 03, Issue 06 


2.7. Simple Algorithm 


Semi-Implicit Method for Pressure Linked Equa- 
tions is an essentially a guess-and-correct procedure 
for the calculation of pressure on the staggered grid 
arrangement. 


STEP 1: Solve discretised momentum 
equations 


Qi; Uz; = LAnyth, + (Pia; — Pi; )AiLy + 42; 


STEP 2: Solve pressure correction equation 


‘ oniis , ‘ , 
Qi jPijg = G1 jPi-vyj + Giza jPiesy + @ij—-1Pij—-1 
+ Qi j41Pije. + b: 


STEP 3: Correct pressure and velocities 
Pig = Piz + Pi; 


Uz = UZ; + a; ;Mi-+.35 — Pid 


STEP 4: Solve all other discretised 
transport equations 


Qi O77 = Aj-1 j;Oi-1 jg F Aja ;Oi41 ; 
+ Q; 5-19; 5-1 


Convergence? 


STOP 


Chart 1: Simple Flow Chart 





TABLE 1. Thermophysical Properties of Fluid 
and Nanoparticles 


Proper- Ag Cu =Al203 TiO, Water 
ties 

Cp 235 383 =6765 686.2 4179 
Ro 10500 8954 3600 4250 997.1 
K 429 400 46 8.954 0.6 
Beta 5.4 16 0.63 2.4 21 


3. Results and Discussion 


In the partially heated heater analysis of buoyancy 
flow induced natural convection by simulation work 
using fundamental Navier Stock equation and tech- 
niques used was of computational fluid dynamic to 
discretised the velocity, temperature, pressure equa- 
tion. The base case was taken with A=1, Ra=10°5, 
g=0.05. 
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Figure 4(a). Cu-Water 
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Figure 4 (6) Ag-Water 
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Figure 4 (c) TiO?- Water 





Figure 4 (d) AhO)-Water 


FIGURE 4. Streamline (Left) and Temperature 
Variation (Right) (a) Cu- Water (b) Ag- Water 
(c) T10, Water (d) AlO.- Water 


TABLE 2. Thermophysical Properties of Fluid 
and Nanoparticles at Ra =10°5. 
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FIGURE 5. Variation of average Nusselt number 
corresponding to different volume fraction (cof) 
using different nanoparticles, at defined Rayleigh 
Number Ra =10°S. 
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FIGURE 6. Variation of average Nusselt number 
corresponding to different volume fraction (cof) 
using different nanoparticles, at defined Rayleigh 
Number Ra =10°3. 


TABLE 4. Thermophysical Properties of Fluid 
and Nanoparticles at A= 0.5, 1 & 2 


C.o.f Ag Cu Al203 TiQs 

0.05 6.4532 6.4802 6.2841 6.1758 
0.1 7.4856 7.5296 7.1629 6.9169 
O.15 8.5856 8.6338 8.1403 7.7219 
0.2 9.7729 9.8157 9.2358 8.6035 


TABLE 3. Thermophysical Properties of Fluid 
and Nanoparticles at Ra=10°3. 


C.o.f Ag Cu Al203 TiOQs 

0.05 6.1713 6.1711 6.1411 6.0192 
0.1 7.1024 7.1019 7.0354 6.7971 
0.15 8.1425 8.1415 8.0302 7.5851 
0.2 9.3118 9.3103 9.1434 8.4837 
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C.o.f A=0.5 A=l A=2 
0.05 5.9847 4.8092 3.0462 
0.1 6.9795 5.5897 4.1686 
0.15 7.9878 6.4045 4.7487 
0.2 9.0391 7.2730 5.1589 
Figure 4. [a-d] shows Streamline (Left) and 


Temperature Variation (Right) where a compari- 
son between Cu-water, Ag-water, TiO:-water and 
Al,O3-water using base conditions is explained. 
A=1, Ra=10°5, 6=0.05, iteration =100, convergence 
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—a— A=0.5 





0.04 0.06 0.08 0.10 0.412 014 O16 0.18 0.20 0.22 
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FIGURE 7. Variation in value of average Nusselt 
number corresponding to different values of @ at 
defined A=0.5, 1 & 2. 


of result =O and mesh define for the rectangu- 
lar cavity as nx=50; ny1l=30; ny2=60; ny3=40; 
ny=ny |+ny2+ny3-2; 

Figure-5. Variation of average Nusselt num- 
ber corresponding to different volume fraction (cof) 
using different nanoparticles, at defined Rayleigh 
Number Ra =10°5. The figure here shows the mono- 
tonically increase in heat transfer corresponding to 
increasing volume fraction for all the mentioned 
nanofluids at three values of Rayleigh number. For 
Ra=10°5 there was negligible changes found for Ag 
and Cu and for further decreases value of Ra the Cu 
and Ag shows almost the same behaviour of increas- 
ing heat transfer and gave same mean Nusselt num- 
ber. The lowest value was seen for TiO, because of 
the dominating of conduction mode of heat transfer. 
Therefore, at high Rayleigh number the convection 
mode of heat transfer is more and as we increased 
the volume fraction of nanofluids, the larger differ- 
ence for mean Nusselt number was gained. 


4. Conclusion 


This study was based completely on simulation 
work by using the software MATLAB to investigate 
the effect of any nanofluid with water as base fluid 
inside the rectangular shaped cavity. The impor- 
tant points drawn from this study was Increasing 
the Rayleigh no. and heater scale enhance the heat 
transfer and fluid flow strength without affecting 
other discussed parameters and Different formation 
of heat and fluid flow was attend at different volume 
fraction of different nanofluids. Heat transfer values, 
Nusselt no. values all depend on different nanoflu- 
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ids. In future study, for the rectangular profile and 
with any nanofluids this study can easily. 


Data Availability Statement 


All data that support the findings of this study are 
included within the article (and any supplementary 
files) 


ORCID iDs 


Ajay Trehan © https://orcid.org/0000-0002-7097- 
1094 


References 


Ahn, Chan-Sol, et al. “Experimental, analytical, 
and computational study of natural convection in 
asymmetrically-heated vertical shafts’. Interna- 
tional Journal of Thermal Sciences 170 (2021): 
107131-107131. 10.1016/j.1jthermalsci. 2021. 
107131. 


Akinshilo, Akinbowale T. “Thermal  Perfor- 
mance Evaluation of MHD Nanofluid Trans- 
port Through a Rotating System Undergoing 
Uniform Injection/Suction with Heat Genera- 
tion”. BioNanoScience 9.3 (2019): 740-748. 
10.1007/s12668-019-00653-9. 


Alizadeh, Rasool, et al. “Artificial intelligence pre- 
diction of natural convection of heat in an oscil- 
lating cavity filled by CuO nanofluid”’. Journal of 
the Taiwan Institute of Chemical Engineers 124 
(2021): 75-90. 10.1016/j.jtice.2021.04.067. 


Azad, Mohammad, Dominic Groulx, and Adam 
Donaldson. “Natural convection onset during 
melting of phase change materials: Part I — Effects 
of the geometry, Stefan number, and degree of 
subcooling”. International Journal of Thermal 
Sciences 170 (2021): 107180—107180. 10.1016/ 
j.jthermalsci.2021.107180. 


Chen, Yanjun, et al. “Numerical simulation and 
analysis of natural convective flow and heat trans- 
fer of nanofluid under electric field’. Interna- 
tional Communications in Heat and Mass Trans- 
fer 120 (2021): 105053-105053. 10. 1016/j. 
icheatmasstransfer.2020.105053. 


Ho, C. J., C. Jang, and Chi-Ming Lai. “Natu- 
ral convection heat transfer in PCM suspensions 
in a square enclosure with bottom heating and 


29 


Trehan and Srivastava 


top cooling”. Thermal Science and Engineering 
Progress 25 (2021): 101037—101037. 10.1016/j. 
tsep.2021.101037. 


Kumar, Pradeep, G. Chanakya, and Naman Bart- 
wal. “Investigations of non-gray/gray radiative 
heat transfer effect on natural convection in tall 
cavities at low operating temperature’. /nterna- 
tional Communications in Heat and Mass Trans- 
fer 125 (2021): 105288-105288. 10. 1016/]. 
icheatmasstransfer.202 1.105288. 


Loenko, Darya S., Aroon Shenoy, and Mikhail A. 
Sheremet. “Effect of time-dependent wall temper- 
ature on natural convection of a non-Newtonian 
fluid in an enclosure’. /nternational Journal of 
Thermal Sciences 166 (2021): 106973—106973. 
10.1016/.14jthermalsci.2021.106973. 


Ma, Yuan, et al. “MHD forced convection of 
MWCNT-—Fe304/water hybrid nanofluid in a par- 
tially heated v-shaped channel using LBM”. Jour- 
nal of Thermal Analysis and Calorimetry 136.4 
(2019): 1723-1735. 10.1007/s10973-018-7788- 
4. 


Mahdi, Mustafa S., et al. “Natural convection 
improvement of PCM melting in partition latent 
heat energy storage: Numerical study with exper- 
imental validation”. /nternational Communica- 
tions in Heat and Mass Transfer 126 (2021): 
105463—105463. 10.1016/).icheatmasstransfer. 
2021.105463. 


Okabe, Takahiro, et al. “Spatio-temporal measure- 
ment of natural convective heat transfer on 
melting process using infrared thermography”. 
International Journal of Heat and Mass Trans- 
fer 181 (2021): 121882-121882. 10. 1016/]. 
ijheatmasstransfer.2021.121882. 


Pranowo and Agung Tri Wijayanta. “Numerical 
solution strategy for natural convection problems 
in a triangular cavity using a direct meshless 
local Petrov-Galerkin method combined with an 
implicit artificial-compressibility model’. Engi- 
neering Analysis with Boundary Elements 126 
(2021): 13-29. 10.1016/}.enganabound. 2021. 
02.006. 


International Research Journal on Advanced Science Hub (IRJASH) 


2021, Vol. 03, Issue 06 


Susantez, Cigdem, Aldélio B. Caldeira, and Bruna 
R. Loiola. “Natural convection effects on TNT 
solidification inside a shaped charge mold”. 
Defence Technology (2021). 10.1016/j}.dt.2021. 
7 O10. 


Taghilou, Mohammad and Ehsan Khavasi. “Ther- 
mal behavior of a PCM filled heat sink: The 
contrast between ambient heat convection and 
heat thermal storage’. Applied Thermal Engineer- 
ing 174 (2020): 115273-115273. 10.1016/j. 
applthermaleng.2020.115273. 


Venkateshwar, K., et al. “Influence of Natural Con- 
vection on Design of Thermal Diodes’. [nterna- 
tional Communications in Heat and Mass Trans- 
fer 127 (2021): 105574—-105574. 10. 1016/j. 
icheatmasstransfer.202 1.105574. 


Yang, Dezhi, et al. “Thermal transport of natural 
convection flow of second grade bio-nanofluid 
in a vertical channel”. Case Studies in Ther- 
mal Engineering 28 (2021): 101377—101377. 10. 
1016/}.csite.2021.101377. 


Yang, Pengying, et al. “Simulating natural con- 
vection with high Rayleigh numbers using 
the Smoothed Particle Hydrodynamics method”. 
International Journal of Heat and Mass Trans- 
fer 166 (2021): 120758—120758. 10. 1016/j. 
ijheatmasstransfer.2020.120758. 


© Shivani Srivastava et 
= al. 2021 Open Access. 

This article is distributed 

under the terms of the Creative Com- 
mons. Attribution 4.0 International License 


(http://creativecommons.org/licenses/by/4.0/), 
which permits unrestricted use, distribution, and 
reproduction in any medium, provided you give 
appropriate credit to the original author(s) and the 
source, provide a link to the Creative Commons 
license, and indicate if changes were made. 


Embargo period: The article has no embargo 
period. 


To cite this Article: Trehan, Ajay, and Shivani Sri- 
vastava. “Solving Natural Convection by Navier- 
Stokes Equation in a Closed Cavity using MAT- 
LAB.” International Research Journal on Advanced 
Science Hub 03.06 (2021): 124-130. http://dx.doi. 
org/10.47392/rjash.2021.151 


130 


